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A model for the response of a poroelastic layered soil to an incident Plffe**^* develop^ in a 
previous paper [Saltier etai, J. Acoust. Soc. Am. 79. 1345-1352 (1986) ] is used to pred.ct the 
implex Lndpressure within theupper poroelastic layer. Thepredictionsbothofph^^^^^ 

and attenuation of the slow wave associated primarily with propagation m the pore flmd are 
compared with measurements made with a specially constructed probe microphone. The 
agreement between theory and experiment is good. The predictions for a layered Poroelastic soil 
model are compared numerically with those of a semi-infinite rigid porous sod model and are 
found to differ only at frequencies higher than 1000 Hz. Analysis of the sensitivity to the 
' theoretically predicted propagation constants in the poroelastic soil to the assumed value for the 
bulk rigiditymodulus of thesoil predicts that over theknownrangeofngidj^mod^^^^^ 
possibletoobtainaswitchoverbetweenfastandslowpropagationmodes.mssw.tchover^^ 
at^e lower end of the possible range of values of the shear modulus. It is suggested that probe 
microphone measurements in air-fiUed soils offer a way of measuring Ao^^^^^^^^y ^^^^J^^ 
deduchig the structural parameters required for appUcaUon of the Biot-StoU model to water- 
saturated sediments. 
PACS numbers: 43.28.Fp 

.MTPoniiCTlON published previously.* The results of the attenuation mea- 

INTRODUCTION , suremcnts using the probe indicate that, in many cases, the 

For many years researchers have assumed tnat tne investigated are not acoustically homogen- 

ground surface is acoustically locaUy reacting and have used measurements also reveal extremely slow speeds 

an empirical relationship for outdoor ground unpedanc« propagation of the acoustic wave in porous soU, which 
derived by Delany and Bazley.' These expressions are based ^^^^ common assumption of local reaction. This 
on the assumption that the ground behaves acoustically J^a that there is negUgible difference between 

modified fluid, and they predict the normal surface miped- subsurface acoustic fields due to the incidence of either 
ance and the propagation constant from thebulk flow res^- .^^^ ^aves as has been predicted elsewhere.' 

tivity only. In other instances the ground has been modeled 

as a semi-infinite rigid porous medium,^ as a variable poros- . couSTIC TRANSFER FUNCTION 
ityririd-porousmediom,'-* or as a hardback rigid porous I. ACOUSTIC THAN5>i-t« ruww 

laver^* These models require additional parameters for the The acoustic transfer function of the ground surface 

detcmiinationofthecharacteristic impedance and the prop- may be defined as the nitio of the acoustic pressure at som^ 

agationconstantintheporoussoawhichareneeded.intum, specified depth below the surface to that above the surface, 

todetermine the surface impedance. In a previous pubUca- Using the poroelastic layer moJldescrrt)edprenously.ti^^ 

tion.' the surface of the ground was modeled as that of a acoustic pressure above the surfiice results ftom the madent 

poroelastic soil layer of known thickness overlying a semi- and reflected waves and may be expressed as 

infinite nonporous elastic substrate. The poroelastic layer P^ = iKyI,iB, -^-B,) cos0, . (D 

permits two eompressional wive types <jll«*J»J ^^^I where B, and 5, are the ampUtudes of the incident and re- 

respectivdy. In this p^. the propagation co^^^t waves. L is the bulk modulus of the air. /, is the 

slowwavedeterminedfrom^eporoelasticsamMn^^ p^tionloniant for the incident wave. is the angle of 

layered models are compared to those of the semi-mfimte pwp»s» — r xi,.~,r^fl,«rfnr«sure below the 

retical predictions. Some of these compari sons have been depth z = dis 

P^ = UMC-m^M)tiLp{ilt cos did) +iA2liiC-miM)exp{iliCOse^) (2) 
-|-M',/,(C-ni,3ncxp( -//jcos^irf) -l-L4i/,(C-miilOcxp( -//icosdjrf) . 
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In this notation, the two compressional wave propagation 
constants are denoted by li and /j. The elastic constants of 
the composite, medium are denoted by H, C, and M. The 
plane wave solution for the boundary conditions, including 
shear waves, requires nine plane waves. However, only four 
of the waves contribute to the pore fluid pressure. These are , 
the transmitted and reflected fast and slow waves in the layer 
itself for which the amplitudes are denoted by Ai,A2,A{, 
and A 2, where the prime indicates reflected waves and the 
d/*s are the associated angles of refraction and reflection. In 
Stoll*s notation the m, 's are the ratios of the relative fluid 
displacements to the frame displacements for the fast and 
slow waves. 

In this notation the acoustic transfer function between 
the surface and depth c/ is a complex quantity expressed as 

TAd)=P^io.id)/P^,,(0). (3) 

By calculating the acoustic transfer function for two differ- 
ent depths below the surface, both the pore fluid attenuation 
and phase speed can be determined for the poroelastic layer. 
For depths of z, and (Z2>z,), the frequency-dependent 
attenuation coefiicient may be expressed as 

aio)) = [\n\TAz2)\-^\TAzi)\] 

X201ogio^/(^2-^i) f (^) 

where o) is the angular frequency. The phase speed was also 
calculated from the acoustic transfer function at two depths 
below the ground surface and may be expressed as 

c(<y ) = {tan- * [ Im{ rj/Re( 7; ) ] }- * U2 - )<y . 

(5) 

For the semi-infinite poroelastic treatment, the acoustic 
pressure below the surface may be expressed more simply, 
since there are no waves returning from a lower boundary. 
The pressure at a depth d below the surface is expressed 
easily from Eq. (2) as 

= l4i/|(C- miM)exp(//i cos 6id) 

+ iA2l2iC - mjAf) exp(//2 cos 02d) . (6) 

The corresponding transfer function, attenuation, and phase 
speed are given by Eqs. (3)-(5) after noting that the pres- 
sure above the surface is givm by Eq. (1). The complex 
transfer function for plane waves incident at angle 0i on a 
semi-infinite rigid porous interface may be written simply as 

r, (rf) = cxp[i/, din^ - sin* tf, ) , (7) 

and n = k^/l{ is the complex refractive index of the rigid 
porous medium, kt, being the complex wavenumber within 
the medium. The only permitted wave motion within the 
rigid porous medium represents a limiting case of the acous- 
tic slow wave within a poroelastic medium of the same flow 
resistivity, porosity, and pore structure. Analytical expres- 
sions for may be found elsewhere.* Figure 1 shows com- 
parison of the magnitude of the complex transfer function as 
predicted by each f the models described above. The param- 
eter values used correspond to those tabulated in Ref. 7 for 
sandy 8oiL 
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FIG. 1 . Magnitude of the complex transfer function predicted for a layered 
poroelastic soil (layer thickness 2 m), for a semi-infinite poroelastic soil, or 
for a semi'infinite rigid porous soil. Results to the accuracy of this graph are 
equal. 



II. MEASUREMENTS OF THE ACOUSTIC WAVE IN THE 
EARTH 

The measurements of sound pressure below the surface 
of the earth were made with a specially designed probe 
shown in Fig. 2. It consisted of an AKG microphone and 
preamplifier housed in a cylindrical brass tube. The inner 
cylinders shown were used to increase the mass of the probe 
and to seal any possible leaks. The nose cone contains small 
holes to allow the element of the microphone to sense the 
acoustic field A rubber membrane was stretched over the 
microphone element to protect it from any particles or mois- 
ture that might penetrate the nose cone when the probe is 
inserted into a predrilled hole. A comparison of the probe 
microphone response to a standard AKG microphone is 
shown in Fig. 3. The frequency response of the probe is es- 
sentially flat up to 3 kHz where it begins to roll off to a new, 
less sensitive, plateau near 7 kHz. The frequency response of 
the probe is essentially unchanged when the probe is pointed 
away from the sound source, the normal geometry for mea- 
surementis below the earth's surface. 

Before the probe was inserted into the ground, a drill 
was used to bore a hole 1 in. in diameter to the desired depth. 
The probe was then forced into the hole and highly viscous 




N se Con 

no. 2. The design of a probe microphone for nwasurements beneath the 
toilrai&oe. 
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FREQUENCY (KHz) 
FIG. 3. Comparison of free-field response of probe and standard AKG mi- 
crophones. AKG microphone (- - 0, probe facing source (— ). probe fac- 
ing away from source ( — ). 

oil treatment was poured around the probe near ground level 
to further seal the probe entry. This procedure worked well, 
giving reproducible results for sandy soil and loam. 

For clays, the soil in the vicinity of the probe would 
quite often crack, providing air paths not previously avail- 
able. This cracking gave results which diflFered from one at- 
tempt to the next. It was generally assumed in such condi- 
tions that the measurement givmg the smallest sound 
pressure at the probe was best. Occasionally the holes in the 
tip of the probe became filled with dirt when the probe was 
inserted. This occurrence could only be detected when the 
probe was removed and examined, so several data runs had 
to be rejected after the fact. In sand it was not necessary to 
drill holes as the probe could be pushed easily into the 
ground. Furthermore, greater reliance could be placed on 
individual measurements with regard to the sound path. 

Measurements for the attenuation and phase speed in 
the frequency range 50-2500 Hz have been made at several 
sites. Initially, a hemispherical dish (radius 0.5 m) was 
formed in the ground of mean flow resistivity 500000 mks 
rayls/m and filled w^th sand which had a measured flow 
resistivity*® of 85 000 mks rayls/m. Among the other sites, 
results are presented here of measurements made at a sand 
quarry and a large, flat area of dredged sand which had been 
weathered over several years and which supported sparse 
growths of lichen and vegetation. Table I lists the phyacal 
properties measured at various locations on this site. The 
acoustic attenuation at all sites was determined in the folio w- 

ing manner. 

A reference acoustic spectrum was recorded usmg 
swept tones transmitted through a small speaker, received 

TABLE L The physical properties measured at various kxations on the 
dredged sand site. 







Mean flow 




Porosity 


resistivity 




(diy) 


(Nsm-*) 


Location 1 


0.3S2 


116000 


Location 2 


0.314 


70 900 


Locations 


0.2S1 


25 000 


Loeatioo4 


0.284 


41S00 
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by a probe microphone at the sand surface, and processed 
through a scanning filter. The probe was pushed then to a 
depth of rf meters and the swept tone repeated. The attenu- 
ation in dB per meter was calculated from the difference in 
the two spectra per unit length. The measurements were re- 
peated with the probe at several depths up to 20 cm as shown 
in Fig. 4. 

Above the dredged sand, white noise was broadcast and 
the acoustic transfer function was measured by recording 
simultaneously the response of a reference AKG micro- 
phone at the surface and probe microphone at a known 
depth below the surface. Using a spectrum analyzer, the dif- 
ference in the phase was determined for the surface and a 
series of probe depths. This difference is directly related to 
the phase speed of the acoustic wave between two depths. 
The phase speed is calculated using the expression 

v{(o) =(oLx/L<l> , 
where Ax is the probe separation for the two depths and 
is the phase change in radians for the two depths. 

Another technique was used to measure the phase speed 
of the acoustic wave in the hemispherical sand dish. A dual- 
channel digital oscilloscope was used to measure the delay 
between the reference AKG and the probe when a pure tone 
was broadcast. The probe was pushed further into the 
ground and the delay measured again. From the difference 
in the delay times and the change in the probe depth, the 
phase speed was obtained direcUy as the ratio of the distance 
and time measurements. Figure 5(a) shows the phase speeds 
measured in the hemispherical sand dish. 

Typically, in homogeneous sands like that contained by 
the hemispherical dish, the attenuation with depth was uni- 
form. However, probe measurements in weathered soils, 
sand, and clays reveal a variation of rate of attenuation with 
depth that appeared not to be the only result of air leaks. The 
attenuation per centimeter increased with depth below many 
surfaces as shown in Fig. 4. But several measurements at one 
location showed that the attenuation over the first 2-cm 
depths was greater than that of any other equal spacing be- 
low those 2 cm. In these instances, a load-bearing sand crust 
of approximately 1 cm thick had formed. Bdow this crust 
the sand was quite loose. 

,To enable comparison of predictions with these mea- 
surments, the various soils were sampled and mean flow 
resistivities and air porosities were deterniincd. The flow re- 
sistivities were measured with a mpdifibd Leonard's appara- 
tus*® and the air porosities were obtained by weighing a 
known volume on site, drying, and weighing again. 

For each of the three sites investigated the measured and 
theoretically predicted attenuations are shown in Figs. 4, 
5(b), and 6. The predicted values were obtained from the 
rigid porous theory. Figures 5(a) and 7 show the compari- 
son of measured and predicted phase speeds of the induced 
acoustic wave at two sites. 

III. DISCUSSION AND CONCLUSIONS 

In Fig. 1. the magnitudes of the acoustic transfer func- 
tions calculated using the semi-infinite and layered poroelas- 
tic models and the rigid frame model differ by less than 1%. 
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The comparison of the acoustic attenuation predicted by the 
three theoretical models described implies that, for most 
cases, the more rigorous poroelastic theory is not necessary 
to predict the acoustic attenuation within outdoor ground 
surfaces. Moreover, the slow wave speeds, both predicted 
and measured in porous soils (Figs. 5 and 7), confirm the 
hypothesis for the explanation of local reaction advanced 
earlier. 

For a realistic range of flow resistivities, the energy asso- 
ciated with the wave propagating primarily in the pores is 
refracted strongly towards the normal to the surface. The 
impedance then becomes independent of the angle of inci- 
dence and regions of the surface, distant from the point of the 
surface at or below which the sound field is probed, do not 
contribute to the r^ected or refracted acoustic fields. This is 
difficult to understand if the wave speed in the earth is taken 
to be the fast or p-wave speed which is of the order of hun* • 
dreds of meters per second, but foUows readily from Snell's 
law if the wave speed in the pores is only tens of meters per 
second as has been measured The locally reacting nature of 
the ground surface is a natural consequence of the porosity 
of the upper layer* The high attenuation of the penetrating 
acoustic wave also explains why the actual seismic profile of 
many grounds may be ignored when considering sound 
{Propagation near their surfaces. 

For typical flow resistivities, surface layers with thick- 
ness greater than 0.2 m will behave acoustically as though 
they are infinitely thick. Furthermore, as has been demon- 
strated in both this paper and the previous one,^ the bulk 
moduli of typical soils are sufficiently high that, acoustically, 
they may be regarded as rigid. These results support the ap- 
plicability of relatively simple models for the acousti<»l 
properties of ground surfaces as have been developed else- 
where.^ When considering acoustically induced soil particle 
motion, however, the poroelastic models are essentiaL 

It is noticeable that the measured attenuations per centi- 
meter ver the uppermost 2 cm in dredged sand differ oon- 
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siderably from those predicted using the measured mean 
flow resistivity and porosity (Fig. 6). On the other hand, 
adjustment of the flow resistivity, increasing its value to fit 
the lowest frequency data point, enables tolerable agreement 
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FIO. 5. (a) Measured aad predicted phase speed in a hemispherical sand 
(h>M#«tiM^aiidpfi^Wffdanfflnatinn in a hemispherical sand di^ 
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FIG. 6. Measured and predicted attenuation in dredged sand. Theoretical 
calculations based on Eq& (4) and (7) and Ref. 2 using the specified values 
of porosity (O) and flow resistivity {Rs)- Measurements made with probe 
at 2 cm for four different locations at the dredged sand site. 



between theory and data over the whole frequency range of 
> measurement Furthermore, use of these adjusted values of 
.flow resistivity produces theoretical predictions of phase 
speed that are in good agreement with.ineasurements (Hg. 
This suggests that the measured flow resistivity is likely 
be in error for such sandy soiL This may be a consequence 
the invasive technique used, by which a sample is re- 
ved, thus loosening ills structure and reducing its appar- 
~:t flow resistivity. 

Given the otherwise good agreement that has been 
shown between measured and predicted attenuations and 
jphase speeds, the possibility arises of using acoustic mea- 
surements of the kind reported here as an in situ technique 
for measuring flow resistivity. This would provide useful 
.^ta both for outdoor sound propagation and for soU 
ttce. 

It should be noted that the assumption has been made 
that the acoustic wave in the pores is the slowest of the two 
^mpressional waves possible in poroelastic soiL Indeed, for 
Che soils studied here this is the case. However, it is not neces- 
sarily always the case. Figure 8 shqws the result of reducing 
i diear modulus while keeping the remaining parameters 
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FIG. 7. Measured and predicted phase speed in dredged sand. Measure- 
ments were made between 2- and 4-cm depths. Predictions made as in Fig. 6. 
Error bar is deduced from attenuation measurements made at these depths. 

constant as listed in Table XL A switchover between the two 
compressional modes is predicted at high frequencies as the 
shear modulus is reduced. Consequently, it would be possi- 
ble in a loosely packed soil with high compressibility and low 
rigidity for the frame wave to be the slowest in some range of 
frequencies. This is the case, in fact, for many polyurethane 
foams, as has been observed elsewhere.*.' 

A final remark about the acoustic wave in air-filled soils 
cpncems the ease with which it may be measured and pre- 
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no. 8. Effect of redttdttg shear modulus to 1.38 x 10* dyn/cm' on fast (— ) 
and slow ( — ) wave speeds in poroelastic soiL 
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TABLE IL ParBmctm used for Fig. 8 calcuUtions. 



Parameter 



Sandy soil 



Flow resistivity (N s m"^) 
Porosity (diy) 

Frame bulk modulus (N/m^) 
Grain shape factor 
Shape factor ratio 

Grain density ( kg/m' ) 



78 300 
0.4 

8.27x10^ 
1 

0.75 
2 650 



dieted at audio frequencies, in comparison with its counter- 
part in water-filled sediments whose existence was not veri- 
fied at ultrasonic frequencies until fairly recently. The fact 
that the acoustic wave dominates the energy transfer across 
the boundary of an air-filled porous material, whose frame is 
relatively rigid, means that measurements of sound reflec- 
tion and transmission at the surface of fibrous absorbents 
have demonstrated its existence since the time of Lord Ray- 
leigh." Consequently the existence and properties of the 
acoustic slow wave in a porous material were known well 
before their prediction by Biot.'* Furthermore, values of 
structural parameters, introduced by the Biot theory, that 
are independent of the saturating fluid, may be determined 
by studies of air-filled rather than water-filled media. 
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